Self-force on a scalar particle in a class of wormhole spacetimes 
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We consider the self-energy and the self-force for scalar massive and massless particles at rest in 
the wormhole space-time. We develop a general approach to obtain the self-force and apply it to 
the two specific profiles of the wormhole throat, namely, with singular and with smooth curvature. 
We found that the self-force changes its sign at the point where nonminimal coupling £ = 1/8 
(for massless case) and it tends to infinity for specific values of £. It may be attractive as well 
repulsive depending on the profile of the throat. For massless particle and minimal coupling case 
the electromagnetic results are recovered. 
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I. INTRODUCTION 

qq , The wormholes are topological bridges which connect different universes or different parts of the same universe. 
These kind of tunnels in space-time have been appeared in different contexts of physics: in the analysis of black 
hole backgrounds [H, 0] , as an idea to construct " charge without charge" and " mass without mass" [!, 0| , and 
as a possibility for time machine p, . The wormholes require an amount of exotic matter which breaks energy 
qh| conditions of a solution of Einstein equations. There are many exotic ideas how to produce exotic matter and 
how much exotic matter one needs to make possible existence of wormhole. Carefull discussion of wormhole's 
^q' geometry and physics may be found in the Visser book Q, in the review by Lobo Q, as well as in paper Q 

i 1 related with astrophysical implementation of wormholes. Among interesting publications we can mention Rcf. 

flpj l which analysed numerically the process of the passage of a radiation pulse through a wormhole and the 
subsequent evolution of the wormhole. It was shown that the wormhole is unstable and it is transformed into 
a spacetime with horizon. The analysis was made for normal as well as exotic matter pulse formed by scalar 
fields. 

. It is well known that a particle in curved spacetimes may interact with the gravitational background by specific 
interaction due to self-force The origin of this self- force is associated with nonlocal structure of the field, 
the source of which the particle is. The self- force may be the unique gravitational interaction on a particle as it 
happens for particles in the space-time of a cosmic string [12[ • In this case there is no gravitational interaction 
ON . particle with cosmic string, but nevertheless there exist the repulsive self-interaction force. In contrast to 
standard self- interaction Dirac-Lorentz force [l3j . the self-interaction force in curved space-times depends on 
all history of the particle and it is usually non-zero even for a particle at rest. For a particle at rest it may be 
found as coincidence limit of the renormalized Green function [lj| . A detailed discussion of the self- force maybe 
found in reviews 

sa' 



In a recent paper [17| the self-force for an electromagnetically charged particle at rest in the static wormhole 
background was analyzed in detail. The general expression for self-energy for arbitrary profile of the throat was 
obtained. It was shown that the particle is attracted by the wormhole and this effect may has astrophysical 
applications. For a specific profile of the throat the result was confirmed by Linet in Ref.[18|] using a different 
approach. There is also another approach for this question which was developed by Krasnikov in recent publica- 
tion for a specific profile of the throat. The difference of results connects with understanding the self- force 
itself. The self-force for scalar particle reveals peculiarities [2(1 Hi], HH due to nonminimal coupling of the 
scalar and the gravitational fields. For example, in Schwarzschild spacetime the self-force on a particle at rest 
is zero for minimal coupling [2(|. In the present paper we analyse in detail the self-force on a scalar particle at 
rest in the background of wormholes. We found that the self-force has crucial dependence on the nonminimal 
coupling constant £: it is zero for £ = 1/8 in massless case and it tends to infinity for specific values of £. 

The organization of this paper is as follows. In Sec. II we develop our approach and consider the origin for 
the divergence of self-energy for specific values of £, from the point of view quantum mechanics. In Sec. Ill we 
consider the self-energy and self-force for massive and massless cases for two specific profiles and for a general 
profile of the throat. Section IV is devoted to the discussion of the results. Throughout this paper we use units 
c=G = l. 
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II. APPROACH 



Firstly let us say some words about the background under consideration. We use the following line element 
of the spherically symmetric wormhole space-time 

ds 2 = -dt 2 + dp 2 + r 2 (p)dn 2 , (1) 

where the profile function r{p) describes the shape of the throat and dVl 2 — d9 2 + sin 2 9 'dip 2 . The variables 
are defined in the following ranges: t,p € R, 9 e [0,7r] and ip £ [0, 2tt]. The radius of the throat is defined by 
a = r(0) and r'(0) = 0. The nonzero components of the Ricci tensor and scalar curvature read: 

2r" 

R p = - — 

V r > 

R e = R % = + r 2+ ^ 

R _ 2(-l + r' 2 + 2rr") 

The three dimensional section corresponding to constant time of this space-time is conformally flat. Indeed, let 
us consider the 3D flat space in spherical coordinates f , 9, p. Thus, we can write 

dl 2 fl = df 2 + f 2 dn 2 . (3) 
Let us choose a new radial coordinate p by the relation r = r{p)e cr ^ p ' > with 

a = ± f P ±--]nr{p). (4) 
J r[x) 



Using this coordinate system we obtain 

1% 



dl 2 = e 2 °{dp 2 + r(p) W) = e 2 °dl 2 wh . (5) 



Therefore g™ k h = e~ 2<7 g{l and the 3D section is conformally flat. For this reason it is expected [24| that the 
self- force is zero for £ = 1/8 and m = as will be shown by manifest calculations. More information about the 
wormhole's space-time may be found in book and review 0] • Some figures for specific profile of throat are 
in Ref. [25[, and the physics in this spacetime is discussed in Ref. [26[. 

Let us consider a massive scalar field, <f>, with scalar source, j, which lives in the wormhole background with 
non-conformal coupling £. The action consists of two part, the first one is for the field itself and the second one 
describes the interaction of the source, a scalar charge e, with the field and is given by 

S = -F- / (^'^ + ^ + ™ 2 2 )V=gd 4 x + f j^V^^x. (6) 



8tt . 

The variation of this action with respect to the metric gives the energy momentum tensor of the field with 
contribution due to the interaction field with charge 



T/.r jog^u + ($;v.4>;u - ^9^-a^ a - \^a 2 g^ 2 ^ + J- (G^(j) 2 + g^D(/) 2 - (j) 2 .^) , (7) 



while the variation with respect to the field gives the equation of motion 

(□ — £R — m 2 )(j) = -iirj, (8) 

with scalar current 

j(x) = ej SW{x-x p {t))-j=. (9) 

We consider only the case in which the particle is at rest at the point x p in the wormhole space-time. This 
means that there is no dependence on the time and the equation of motion for the field has the following form 

(A - to 2 - $R)<t>(x) = -47Tj = -^=5 {3) (x - x p ). (10) 
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From the general point of view the full energy of particle reads 

E = - ( T^dY?, (11) 



where E is 3-surface with £^ as a normal. The spacetime under consideration possesses the time-like Killing 
vector = Sq . Thus choosing the hypersurface of constant time we obtain 

E= - [T m ^d 3 x, (12) 



where for the static case under consideration we have 

Too = -M + -L {4>-,iV + (£R + m 2 )0 2 ) - J- + 0A0) . (13) 

Integrating by part and taking into account the equation of motion (|10[) we get 

i r . , , — i 



E = - I J^V=9d 3 x=- J J j(x)G(x 1 x') 3 {x l )^-g(x)^-g{x')d 3 xd 3 x' 
e 2 

= — G(x p , x p ). 

The Green's function G(x,x') which we need for calculation self-force is divergent in the coincidence limit, 
x 1 — > x. Several methods to obtain a finite result are known. The most simple way is to consider total mass 
as a sum of observed finite mass and an infinite electromagnetic contribution. Usually this procedure is called 
"classical renormalization" because there is no Planck constant in the divergent term. Dirac [13j suggested 
to consider radiative Green's function to calculate the self-force. Since the radiative Green's function is the 
difference between retarded and advanced Green's functions, singular contribution cancels out and we obtain 
a finite result. There is also an axiomatic approach suggested by Quinn and Wald (2?j|. In the framework of 
this approach one obtains finite expression by using a "comparison" axiom. This approach was used in Refs. 
[H, H3| for a specific space-time. 

There is a problem with renormalization for massive uncharged particle (see Ref. [30|). The self- force for 
charged particle is obtained by classical renormalization of particle's mass, m — > m + Ae 2 , where A is some 
infinite constant. Therefore for massive and uncharged particle we have to renormalize mass by using the same 
mass of particle. To solve this problem in the Ref. [30] has suggested to consider particle as Schwarzschild 
solution in external gravitational field. In this way it is possible to obtain finite result in lower power of mass. 
More rigorous derivation has been done in Ref. [3l| . 

We will use a general approach to renormalization in curved space-time [32| , which means subtraction of the 
first terms from DeWitt-Schwinger asymptotic expansion of the Green's function. In general there are two kinds 
of divergences in this expansion, namely, pole and logarithmic ones [33| . In three-dimensional case which we 
are interested in there is only pole divergence, while the logarithmic term is absent. The singular part of the 
Green's function, which must be subtracted, has the following form (in 3D case) 

1 (A 1 / 2 



4tt [ V2ct 

where a is half of the square of geodesic distance and A is the DeWitt-Morette determinant. If we take 
the coincidence limit for angular variables then these quantities are easily calculated using the metric ((T|): 
(7 = (p — p') 2 /2 and A = 1. Thus to carry out renormalization we should subtract the singular part of the 
Green's function, which has the following form: 

Qsing = J__L^ + ™_ 
4n\p-p'\ An 

This approach was used many times in different curved backgrounds (e.g. see [16]). Therefore the self-energy 
U for particle at rest at the point x p has the following form 

2 2 

U(x p ) = E - 7 rG rea {x pi x p )= e - lim {G(x p ,x')-G sm9 (x p ,x')}. (14) 

For simplification of notations we will denote hereinafter the position of the particle by x — x p . 

In order to find the self-energy we have to calculate the 3-D Green function which obeys the equation 

(A-m 2 -gfl)G(x;x') = - (15) 
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and then adopt the renormalization procedure. Due to spherical symmetry we can expand the Green function 
using spherical functions as 

OC ~\~ 1 OO 

G(x;x') = £ £ Y^imA^g^P') = E ^^(cos 7 k(A//)- (16) 

;=0m=-2 i=0 

Here P; are the Legendre polynomials and cos 7 = cos 8 cos + sin 8 sin 0' cos(<^ — ip'). 
The radial Green function,^ (p, p'), obeys the following equation 

Differently from the electromagnetic case [l7| we have a contribution which arises from the non-minimal 
coupling, even for massless case. 

Now, we adopt the approach developed in Ref. [l7| with some modifications. We represent the Green function 
in the following form 

gi = 8{p - pO* 2 (p)*i(p') + 6{p' - p)¥ 3 0/)*i(p), (18) 
where ^1 and \&2 are independent solutions of the corresponding homogeneous equation 

2r' / 2 1(1 + 1) 

with boundary conditions 



+ — - m 2 + v „ ; + £R * = 0, (19) 



and Wronskian condition 



lim * 2 = 0, lim *i ^ 0, (20) 

p — >-\-oo p— >+oo 



^! 1 %-9 2 % = ~. (21) 



It is worthily calling attention to the fact that if we change the function 'J = $/r in Eq. (fl"§)) , we obtain 

*" + (-™ 2 - - t R - y) $ = °- ( 22 ) 

From the quantum mechanical point of view, Eq. (|22p describes a quantum particle in the potential 

r" 

V = £P + — (23) 
r 

Let us consider this potential far from the wormhole's throat, p ^> a. The behavior of the potential crucially 
depends on the profile function r(p). We divide the profiles of the throat in two large classes depending on its 
behavior at infinity. The first class may be called as wormhole without parameter of the throat's length. In this 
case the expansion of the profile function over p — > 00 has a polynomial form 

00 

r = p+J2 b ^\ (24) 

which starts from p~ n with some integer n > 0. The potential far from the wormhole's throat has the following 
expansion 

where £ c = (n + l)/4n. The point £ c is crucial because at this point the potential changes its sign. The smaller 
n, the greater the critical value of £ and vice-verse. The case n = in some sense corresponds to the following 
"weak" dependence of the profile of the throat on p at infinity 

r{p) = p + a + WpT (26) 
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with n > 0. Indeed, in this case there is no critical value of £. This kind of expansion corresponds to n — > in 
the above case and £ c — + oo. 

Another kind of wormholes have a dimensional parameter, t, which describes the throat's length. In this case 
the expansion for p — > oo has the following form 

r (p) = P + a + Cnp n e-T, n > 0, (27) 

and the space-time becomes flat exponentially. The expansion of the potential starts from the following term 

£R+ — = cJ—^p^eS + ■■■. (28) 

Therefore the critical value of £ is 1/4 and it does not depend on n. Therefore, we claim the following statement: 
The critical value for the first kind of throat is £ c = (n + l)/4n if r — p ~ p~ n , and for second kind is £ c = 1/4. 
Obviously the delta function is a short range potential and it belongs to the second kind of wormhole as we will 
see later. In the mixed case the main role is played by the polynomial part of the expansion. 

Let us consider specific examples. For r — \J p 2 + a 2 we have the first kind of throat and n = 1, b% = 1/2, 
the critical value of £ c is 1/2. For the profiles 

r = p coth — + a — r, 
r 

r = ptanh — + a, 

T 

we have correspondingly n = 1,6 = 1/2 and n = 1,6 = —1/2 and therefore the critical value £ c = 1/4. The 
case of a singular potential is a limiting case of shortest length of the throat and it belongs to the second case. 
Therefore, we can expect peculiarities for £ w £ c . Below we will see them in manifest forms. 

Let us consider the case of singular scalar curvature separately. The point is that for the throat profile 
r = \p\ + a, the scalar curvature reads 



R=—S(p) (29) 



a 



and we have the following radial equation 



vi/" + _ (m 2 + l l±ll - ^S(p) ) * = 0. (30) 



r \ a 

Integrating this equation around p = wc obtain the following matching conditions at the throat 

¥(+0)-¥(-0) = 0, 

*'(+0) - *'(-0) = - — *(+0). (31) 
a 

Let us now represent the solutions \& as linear combinations of two independent solutions in each domain of the 
wormhole spacetime, "+" and which correspond to the signs of p. In each domain we find two independent 
solutions 4>± and <f>± with the condition that falls down for p — * oo and the Wronskian condition 

W(4>1,4? ± ) = ^-, (32) 
where A± are constants. Therefore, we have in general 

f a \dp\ + p>0 

1 ~ \ oL<\> l _ +0Ltt?_, p < ' 

_ f aUX + pU* P >{) 

2 ~ \ a 2 _^_+f3 2 ^, p<0 ■ 

The Wronskian condition implies the following constraints on the coefficients: 

1 fl 2 fl i 2 1 
a±P±-P±a± = 

Taking into account the boundary conditions PT|) . then the solutions ^\ and ^2 turns into the following 
forms 



6 



* 2 = al^+^J 2 , 



where 



0i ,p<o, 

^- ,P<0, 

4>i , p > o, 



^+1 wX0Y^J 



f </>+ >P>0, 

In order to satisfy the boundary conditions for 'J, namely linip^oo vf^ = 0, wc have to take a 2 + = 0. For the 
second solution there is no additional condition to leave only one constant. In the space-time without wormholc 
we have the additional condition at the point p = and therefore, the functions must be finite. Here we have 
no origin, there is a bridge starting from some distance of throat. For this reason we have to consider both 
possibilities. Let us consider the specific solution for a second solution for a\_ = 0. It means that we consider 
the function which is symmetric to ^27 an d tends to zero in the mirror spacctimc for p — > —00. These solutions 
have the following form: 



* 2 = I3 2 J 2 - 



(33a) 
(33b) 



In what follows we will consider the symmetric profile of the throat r(—p) = r{p). Taking into account the 
above relations we obtain the radial Green's function in the following form 
1. p > p' > 



2. < p < p' 

3. p < p' and p' > 0, p < 

4. p > p' and p' < 0, p > 

5. pf < p < 

6. p < p' < 



A + w + {<t>\A\) + % i<t>W + 

g\ 2 \p,p')=g\ 1 \p',p) 
W{4>\,4>\) 



^W)= 5l (3) (pV) 

gl B \p,f/) = gl 1 \-P,-P J ) 



WW+{-p) 



g?\p,p')=g ( f\p',P) 

in man 

The second kind of solutions has the following form 



In fact, we have to write out only g^ and gj® in manifest form. 



* 2 - /?^ 2 . 



(34a) 



(34b) 



(34c) 



(34d) 



(34e) 



(34f) 



(35a) 
(35b) 
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In this case the Green functions are given by 
1. p > pf > 



ffpW) = -^+i/'><.';(/<') + 



2. < p < p' 



3. p < p' and p' > 0, p < 

4. p > p' and p' < 0, p > 

5. pf < p < 



^(p')^(-p) 



9? ) {p,t*) = i-<p + {-it)<i>\{-p)- 



1 w 4 



6. p < p' < 



gi%,(/)=gl 0) W, P ) 



4>+{-it)M-p) 



(36a) 



(36b) 



(36c) 



(36d) 



(36e) 



(36f) 



For smooth background we have to set ^ = in the above formulas except for the differential equation for the 
radial Green function. 

From relation (|36e[) we observe that the second solution gives a divergent Green function in the domain p < 

(5) 

because the function g\ contains multiplication of functions which tends to infinity for great p. For this reason 
we have to throw away this solution and consider them as unphysical and thus, the Green function is uniquely 
defined. 



III. SELF-ENERGY AND SELF-FORCE 

Let us firstly call attention to the fact that for scalar particles it is possible to solve the problem concerning 
the determination of the self-force and self-energy in a closed form, only for an infinitely short throat and for 
the profile y/ p 2 + a 2 , in massless case. For arbitrary profile there is a problem to solve this problem for the zero 
mode I = 0. Otherwise, for the electromagnetic case, this problem can be solved in a closed form for arbitrary 
profile. 



A. Profile r = \p\ + a 

Let us first of all consider the simplest profile of the throat given by r = \p\ + a. 

1. Massless case 



In this case we may use the same solutions as the one in Ref. [17[ . We have the equation 

xll , 2r' , 1(1+1) 



which has the following solutions 



-4> = o, 

(a±p) L , 

a 2l+1 (a±p)- 1 - 1 . 
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Taking into account the above formulas we obtain the following expression for the Green function 

1 jft r l + l W 2(1 + 1) - 8£' 



r i+i r //+i 2(/ + l) -8£' 

We use now the first expression above in the Eq. (|16p with the coincidence limit for the angular variables, 
namely, 7 = 0, and thus we have 

G(p; p) = ^2 —T—9i(p> p'), 

2=0 

where p' is the observation point and p is the position of the particle. It is very easy to make the summation 
over I. Doing this, we arrive at the following results 

(P'P^ 4tt p — p' 8irrr' i^rr' ' ' ^) ^ ^ 

ae 2 (l-80 fa 2 \ 
U (P) = l^^M-^J (38) 

for Green function and self-energy, respectively. The first term in the Green function is the DeWitt-Schwinger 
expansion and we subtracted this term in accordance with the discussion in the end of last section. The definition 
and properties of function <f>, 
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^M-4£U£(l-4£ + n)-^) , (39) 



may be found in Ref. [34J. 

The limiting cases p — > oo and p — > gives us the following results 

lim U = - a 4^, 
p^oo 4p 2 1 — 4£ 

]i mU = e2(1 ~ 8$) (ln^ +7 + ^(1 -4£)). 
p^o 4a a 

For £ = 1/8 we obtain zero as should be the case due to conformal flatness of the 3D section of constant 
time. We observe that according with the discussion above the divergence for f = 1/4 appears. Furthermore 
we observe appearance infinite number poles at points £ n = (n + l)/4. The appearance of problems with the 
delta-like potential was noted in literature (35l . l36l . [3?l l38| . In the limit of minimal coupling we recover the 
result for the electromagnetic field given by 

U-£m(i-$). (40) 

Indeed, in the limit £ — > and for m = the equation for scalar field (fTU)) coincides with the equation 
for component A of vector field and the formula for calculation the self-energy (fT4"]) coincides with that for 
electromagnetic field (see Ref. [I?} )• Therefore, considering the e as an electric charge we obtain that in this 
limit we have to recover the electromagnetic case. 

For the second kind solution given by (|36ep we obtain the following expression 

m^F - («) 

We observe that the series over I is divergent because r > a and r' > a. Therefore, as noted above we have 
throw away this solution due to the fact that it is unphysical. 



2. Massive case 



In this case the radial equation turns to 

2r' ( 1(1 + 1) , 2 \ 
<p H <p - t. h TO = 
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and has two independent solutions, given in terms of the modified spherical Bessel functions, I v and K v , as 

1 I T 



where x — m(a + p) and v = / + 1/2. Using these solutions we obtain the Green function 

K v (mr)K v (mr') 



(1) = K v (mr)l v {mr') ma{I v K' v + I'M + (8g - l)I v K„ 
9l V^ 7 2maK v K' v + (8£ - 1)K* 

„(3) _ 1 



2maK v K' v + (8£ - ljif* 



K v {mr)K v {mr') 



(42a) 
(42b) 



where r = |p| + a, r' = |p'| + a. 

Using the addition theorem for Bessel function [3§| we obtain the following formula (r > r') 



1 00 1 

= V(2Z + l)/„(mr')^(mr) = e -™fr-r ) ; (43 ) 



and consequently the first term in <;> represents the standard Yukawa contribution. The second term can 
not be represented in a closed form and we have to calculate it numerically. After renormalization we get the 
expression (p > 0) 



r i , ±^ ma(I v K' v + I'M + (8g - \)I V K„ 

^ren(P,P) - ^ 2, » 2 maK v K' v + (8f - l)Kl 



Kl{mT \ (44a) 



U{p) = 27re 2 G re „(p,p). (44b) 

At the beginning we discussed the divergence at £ = 1/4. Let us consider in manifest form the first term 
(I = 0) in the renormalized radial Green function 



. ) ( II 



fJi 

It has the following forms 



ma{I v K' v + I'M + (»C - 1)^^ 



2maK v K' v + (8£ - 1)^ 



_ e -2m(r-a) (3 _ g£) + g-2mr (3 _ g£ + 4fflTO ) 

ffo ' re ™ ~ 4mr 2 (l + am~4£) ' ( ' 



2r 2 l-4£' 



m->0 



(46) 



for massive and massless cases, respectively. We observe that there is no singularity for massive case at point 
£ = 1/4. Indeed, for £ = 1/4 we have 

(1) _ _ e -2m(r-a) + g-Smr^j + 4q7tj) 

3o,re„ - 4m 2 ar 2 ' 

But this expression blows up for massless case because the expansion gives us 



-,(!) 



1 



2mr 2 



Therefore in the massive case the expression is no longer singular at point £ = 1/4. The singularity appears at 
point 1/4 + ma/4. The next term with I = 1 will show a singularity at point 1/2 + (ma) 2 /4(l + ma) and so on. 
In general we observe that divergences appear as solutions of the following equation 

2xK' v (x) + (8£ - l)K u {x) = 0. (47) 

For massless case (x = 0) there is general solution of this equation, £„ = (n + l)/4, where v = n + 1/2, and we 
obtain infinite number of poles. For massive case we have infinite number of solutions, too, but it is impossible 
obtain solution in closed form for arbitrary n. From above consideration we may conclude that nonzero mass 
of field gives correction for obtained for massless case. 
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Let us consider the convergence of the series for the Green function, that is, we have to consider the expressions 
for v — > oo and fixed r. With this aim we use the uniform expansion for the Bessel function [40j which is valid 
for great index. We suppose that mr/u and ma/v are constants and use the uniform expansion for functions 

I v (yz) and K v {yz). We make expansion for vgfl m over v~ x — > and then make the expansion of the obtained 
expression over v — ► oo because each term of the expansion depends on v through mr jv and ma/v. Doing this, 
we obtain the following result 

- iJ f H + I 2 ™ 2 (" 2 - » 3 > + <1 - 15* I 4 " 2 ™ 2 

+ 2m 2 (a 4 m 2 - r 2 - m 2 r 4 - a 2 (l + 2m 2 r 2 )) ( + 2m 2 (r 2 -a 2 ) C 2 + C 3 ] + ■■■}, (48) 

where £ = 1 — 8£. We observe that for r > a the series is always convergent. For r = a the series is still 
convergent only for £ = 1/8. Therefore, the energy at the throat, r = a, is divergent for any case, except for 
£ = 1/8. The numerical simulations of the self-energy are shown in Fig. [1] We note that the massive field will 




FIG. 1: The numerical simulation of the self-energy of a massive scalar field for £ = and for different parameters 
ma — (thick line), ma = 0.1 (middle thickness) and ma = 1 (thin line). In the figure at right we show the numerical 
simulation for ma — 1 and for different parameters £ = 1/10 (thick line) up to £ = 1/3 (thin line). For £ = 1/2 it tends 
to infinity as should be the case. 

produce a self-force which is localized close to the throat. It falls down exponentially fast as e~ mr far from the 
throat. This behavior is in agreement with Linet result |4l| . 



B. Profile r = \J p 2 + a 2 

Because the mass of a field leads to suppression of the self-force by factor e~ mr we consider the massless case, 
7tl = 0. The radial equation is written as 



2p ( 1(1 + 1) 2£a 2 



= 0, 



( j o 2 + a 2 ) r \p 2 + a 2 (p 2 + a 2 ) 2 / 
which has two linearly independent solutions (// = v2?)i namely, 4>\ and (jr+, given by 

0+ = c+P-^(z), (49a) 

4>\ = c+Qf(z), (49b) 

with Wronskian 

Wi^^l) - ic+c+^e^. (50) 

Here Pf and are the Legendre functions of the first and second kinds, and z = ip/a. Taking into account 
the above relations we obtain the following expression for the Green function 
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As expected, for £ = 1/2, that is for = y/2£ — 1, we have a divergence as T(l — /i) in the second term of 
the sum in (fSTj) for Z = 0. Indeed, for Z = and /i = lwe have 



P^z) = (52) 
V-z + 1 

Q5(») = — 1 (53) 
— ly z + 1 



and the first term in the scries reads 



- iP -\z')Ql(z) - l^-A Q l^)Ql( z ). (54) 

Therefore, the Legendre's functions are finite but the gamma function is divergent at p — 1. We noted this 
critical value of £ from the quantum mechanics point of view. Let us consider this moment in manifest form 
the calculations starting from the beginning. For Z = and /j = lwe obtain two solutions 

Ml _ X ^2 _ I (55) 



i1 $ ,o 1 



(56) 



where x — p/a. It is easy to see that the matching conditions give the following solutions 

r + ^±^, p> o 

The first solution falls down for p — ► +oo and the second solution falls down for p — ► — oo. Therefore, we will 
get the following expressions 

*i = 

* 2 = ^ 



x~ 



and the Wronskian is obviously zero. Because the solutions have the Wronskian in the denominator we obtain 
a divergent result. 

Unfortunately, there is no closed expression for the above series for arbitrary £. Nevertheless, for the specific 
case £ = 1/8, we may obtain the Green function and the self- force in manifest form. Indeed, in this case p, = 1/2 
and Legendre functions are expressed in terms of simple functions [34j , which helps us to obtain the following 
expression for the Green function 

G(P ' p,) = Ana (p - p')(l + 2^V 4 (1 + x' 2 )i/4 ' (57) 



where p = x + \J x 1 + 1 and x = p/a. Then, we renormalize the Green function and take the coincidence limit 
and obtain: 

G-> ) ,) = [ G -^ 7y ]=0. (58) 

Therefore, as expected, the self- force for £ = 1/8 for this profile of throat is zero. We will confirm this result by 
numerical calculation below. 

To obtain the expression for arbitrary £ we perform WKB analysis of the radial equation 

<t>" + 7V - (^7^ + <t> = ( 59 ) 

and represent the solution of this equation in the form 

<p = e s . (60) 
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Thus, we obtain the following equation for S: 



where v = I + 1/2. The next step is to expand 5 in the following power series in ^ : 



(61) 



J2 »~ nS n- 



(62) 



n=-l 



As noted in Ref. [TtJ we have to consider the term with I — separately. In this case the equation ([55)) 
simplifies considerably (<p stands here for the zero mode only): 



2r' 



ip' - £Rip = 0. 



The general solution of the above equation for our specific profile reads 

cp = C\ cos(p arctan — ) + C2 sin(p arctan — ). 

a a 

Thus, we have the following solutions 



(63) 



Vl = 



2tan^ 



( ki cos(/Lt arctan — ) + ki sin(/i arctan — ) 
\ a a 

7r r p nu . . p } 

— < cos (a arctan — — cot — smu arctan — ) > 



2sin^ 



^ sin/z 



dp 



(64) 



with Wronskian 



W(^>¥>+) = ~4( fc i +fc 2 tan^). 

2 



Therefore we obtain expression 



So 



m cos(2/z arctan I ) — cos(7T/Lt) 



2a/i sin np 



(65) 



(66) 



which does not depend on ki and k^. 

Now we consider terms with I > 0. Substitution of (f6"2"| into [) yields the set of expressions for the functions 
S n . General solution of the first four equations of this chain reads 

= ±-, 

r 



r' 1 



s' 2 = 

^ = 

Sa = 



= ± 



C 

"4 
8 
16 



« r 1 

r H 

2r 2r 
rr (3) + 2r ' r " 



r< 4 V 2 + 2r" 2 r + 4rV 3 V + 2r'V 
2rV" + rr (3) 



128r 



[r' 2 + 2rr"-l] : 



g[r' 2 + 2rr"-l 



-2rV 3 - 10rr( 3 V 2 - r( 8r" 2 + 7rA 4 A r' - r 



(8r" 2 + 7rrW) r' - r 2 (8rV 3 > + 7V 5 )) 



[r' 4 + 2 (2rr" - 1) r' 2 + 4rr" (rr" - 1)] ' 



+ 1 

128 

_ £ 
32 

We observe that (i) all S' k with k > 1 are proportional to C = 1 — 8£, and (ii) all S^j. are full derivative as in 
the case of £ = 0. 



2rV 2 + Arr^W + r ( 2r"' z + rr 



r.(3) r i 



.,"2 



,(4) 
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Now we are in position to calculate the Green function. First of all we calculate the Wronskian and found 



A + = W{4>\AX)r 



1 A? \J2 _ c S +2 +S* 



(S'+ 2 - S' +1 )r 2 



By using the formulas above we obtain 

s« (1 W) = 



Then, we change summation over n and I and get the following result 



(67) 



(68) 



°° 2v 1 1 °° 1 



where 



/o(6) 



+ o(p-p'), 



= -ii n qV-g + o( p - p0) 

fk(b) 1. 3 , 

7 ; w „ = -Ch(k, -) + 0(p-p). 

VKpHpO r 2 

The functions jfe may be found by using simple code in package " Mathematica" . There is no general form of 
these coefficients for arbitrary index but for numerical calculations we need only for some first ones (see Ref. 
[l7T|). The first four coefficients have the following form in coincidence limit 



k{p,p') = 1, 
h(p,p') = -C 



p -1 + r' 2 + 2rr" , -1 + r' 2 + 2rr" 



J2(P,P) = ~C 

3( 2 



8r 



-l + r' 2 + 2rr" 



■dp = -C- 



8r 



(p-p') + 0((p-p') 2 ), 



M/>»P) 

Therefore we obtain 



2/y 



ir" 2 + rrW 



^2- (r' 2 + 2rr" - l) 2 - ^ ( 2r"r' 2 + 4rr< 3 V + r f 2, 
128 v 16 V V 

cos(2fi arctan £ ) — cos(7r/x) 



1 

47T 



1 j-l + lY,^(2k,hj 2k (p,p) 



p — p r r 



fc=l 



2a/i sin 7r/x 



where each term j 2 k is proportional to £ = 1 — 8£. Here (s,p) is t ne Hurwitz zeta function (see, for example 
(3^|). After regularization we arrive at the following formula for the self-energy (p = \/2£) 



U(p) 



1 , 1 Y^V /07 V / \ 1 cos(2^i arctan f) - cos(tt/z) 
- + - VCi?(2A;,-)j2fc(p,p) + 



fc=l 



2a/i sin7T|U 



As expected it is zero for £ = 1/8 and it is divergent for £ = 1/2. Far from the throat we obtain 



U 



ap 



2p 2 tanirp 



(69) 



(70) 



By numerical analysis it is enough to take into account only two terms of the series above. In fact we use half 
of sum the first two terms. The numerical simulations are reproduced in Fig. [2] 
For arbitrary profile of the wormhole we have the following formula 



U(p) 



\+ l -T.^{2k i -) 32k {p,p)+g ( i\p) 
fe=i 



(71) 
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FIG. 2: The numerical simulation of the self- force on a massless scalar field for profile r = \J p 2 + a 2 for different 
parameters from £ = (thick line) up to £ = ^ (thin line) . For £ = | it is zero and for £ = | it tends to infinity. 



with the same j 2 fe a s above and 

fl£W) = -^+(^+(^ + ^-(4 + 4) ^W^(P), (72) 

where A + = W + {tp\,ip\)r 2 {p). The functions tp^ 2 are the solutions of the equation 

2r' 

ip" + — iff - £Rip = 0. (73) 
r 

Unfortunately, differently from the electromagnetic field case, there is no general solution of this equation for 
arbitrary £ and r. For £ = 1/8 it is easy to find a general solution of this equation by using the conformal 
flatness of the equation. They read 

1 1 fP dy 1 _ rp dy 

<p L = —=e J , ip z = —=e J , (74) 

with Wronskian W^ 1 , ip 2 ) = — 1/r 2 . For £ ^ 1/8 we may only make conclusion about behavior of the self-force 
far from the wormhole's throat. Indeed, changing function by the relation ip = w/r we obtain 

w" - + £Rj w = 0. (75) 

Therefore, for great distance we have simple equation w" = 0, with two solutions 

wi = C2/Q, W2 = ci. (76) 

The Wronskian corresponding solutions is 

W{<p u <P2) = -^. (77) 
It is not difficult to show that the solutions with next corrections are 

PI = * (l + ^- + 0{p- 2 ) 
VI = 

and therefore we obtain for great p 



v , = _l(l + 0(p- 2 )), (78) 



V V 2c l </?2/0 

But we can not make any conclusion about sign of these expression. 
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IV. DISCUSSION AND CONCLUSION 



In this paper we considered in details the self-interaction on a scalar particle at rest in the wormhole space- 
time with non-minimal coupling with curvature. The main peculiarities of the self-force on a scalar field are (i) 
mass of field and (ii) nonminimal coupling £. We consider a particle at rest and for this reason all equations 
become effectively three dimensional, because they touch only spacial part of wormhole space-time which is 
conformally flat. For this reason for £ = 1/8 and massless field we expect that the self- force is zero [241 ]. Our 
calculations confirm this result, the self- force zero indeed in all considered above special examples. For £ < 1/8 
the scalar particle is attracted to the wormhole and for £ > 1/8 the particle is repelled by the wormhole. For 
£ = 1/8 the self- force is zero. In the space-time of a black hole [20] one has different behavior of the self-energy, 
in which case it is proportional to £ and the self-force is zero for minimal coupling £ = 0. 

The self-force for scalar massless particles reveals peculiarity for specific values of the nonminimal coupling 
£ c . The energy has a simple pole, (£ — £ c ) _1 , at this point. Because the self-energy is defined in terms of the 
three dimensional Green function we may analyse this pole by using analogy with the scattering theory [42| . 
The combination V — £i? 4- r" jr plays the role of a potential for the wave function in non-relativistic quantum 
mechanics and all information about boundary or scattering states is encoded in the Green function. This 
potential tends to a constant at the wormhole's throat, p — 0, and it falls down to zero far from the throat. The 
critical value depends on the shape of the throat. If the space-time far from the throat differs from Minkowski 
spacetime as b n p~ n , then the critical value £ c = (n + I) /An and the potential, V = nb n {^ c — S,)p n 3 changes 
its sign in this point. If the spacetime goes to Minkowski spacetime exponentially fast as c n p n e~ p ^ T then the 
critical value £ c = 1/4 and the potential, V = 4c„(£ c — £)r _2 p n ~ 1 e~ p / T , also changes its sign at this point. This 
kind of wormholes possesses a dimensional parameter, r, which may be regarded as length of the throat. The 
profile r = \p\ + a gives singular, delta like potential consentrated at throat which has no longer possess the 
length of the throat. This kind of profile belongs to the second type of throat due to localization potential close 
to the throat. We obtain the same critical value 1/4 in this case by manifest calculations (see Eq. ([38])). 

For wormhole with profile r(p) = y/ p 2 + a 2 the Green function for zero mode, I = 0, is expressed in terms of 
two solutions of radial equation 

„_ o a (l-/x 2 ) 
(p 2 + a 2 ) 2 

with p 2 = 2£. This equation describes the particle in the potential V = — a 2 (l — p 2 )/(p 2 + a 2 ) 2 in one dimen- 
sion. It is well-known [42J that the Green function has the poles for energy of boundary states. The point 
p = 1(£ = 1/2) is critical because the potential changes its sign and boundary states appears or disappears. 

The mass of the field gives additional factor e~ mr and leads to localization of the self-force close to the 
wormhole's throat inside a sphere with the Compton wavelength radius mT x (see Fig. [IJ. The self-force reveals 
singularity too but the point depends on the mass of the field. For example, for a simple profile of the throat, 
r = \p\ + a, the singularity appears at the point £ c = 1/4 + ma/4 (|4"5"1) , where a is radius of the throat. 

We developed a procedure and found general expression (fTTj) for the self- force for general profile of the throat. 

But differently from the electromagnetic field case there is no general solution for zero mode gjj (|T2[) in terms of 
profile function 7'. We would like to note that this relation contains two independent solutions of homogeneous 
radial equation for zero mode at observation point as well as at the origin, for p = 0. Because the fall down 
function ip 2 and its derivative appears at the denominator we have to use irregular solution for this function 
in terms of the scattering theory (see, for example [IH). We observe that zero mode gives main contribution 
to self-force and it depends on the global structure of space-time which is in agreement with consideration the 



self- force in black hole background [4J,[45|]. For example, for profile function r = \J p 2 + a 2 we observe that zero 
mode solutions given by Eq. (|64[) contains the integral over profile function and therefore is defined by global 
structure of the space-time. The numerical evaluations (see Fig. \2§ for this profile show that the self force as 
expected changes its sign at the point £ = 1/8: for £ < 1/8 it is attractive and for £ < 1/8 it is repulsive. 
For £ — > 1/2 the self-force reveals singularity as a simple pole (£ — 1/2) -1 and therefore tends to infinity. The 
self- force has extrema for p a and it is zero at origin. Far from the string the self-energy is given by Eq. I|70p . 
Therefore we may say the same conclusion, the scalar particle will be concentrated at the throat for £ < 1/8 as 
for electromagnetic field case fl7j . 
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